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I. INTRODUCTION 

The problem of a magnetic impurity in a superconductor has been extensively studied, 
but is not completely solved because of the difficulty of treating the dynamical correlations of 
the coupled impurity-conduction electron system together with pair correlations. Generally, 
the behavior of the system can be characterized by the ratio of the Kondo energy scale in 
the normal metal to the superconducting transition temperature, Tk/Tc. For Tk/T^ <^ 1, 
conduction electrons scatter from classical spins and physics in this regime can be described 
by the Abrikosov-Gor'kov theory. [T|] In the opposite limit, Tk/Tc ^ 1, the impurity spin 
is screened and conduction electrons undergo only potential scattering. In this regime s- 
wave superconductors are largely unaffected by the presence of Kondo impurities due to 
Anderson's theorem. p| Superconductors with an anisotropic order parameter, e.g. p-wave, 
d-wave etc., are strongly affected, however and the potential scattering is pair-breaking. 
The effect of pair breaking is maximal in s-wave superconductors in the intermediate region, 
Tk ~ Tc, while in the anisotropic case it is largest for Tk/Tc ^ oo. 

Regardless of the type of the pairing mechanism, the conduction band density of states 
is depleted in the neighborhood of the Fermi surface, resulting in the weakening of the 
coupling to the impurity. On the other hand, increasing the impurity concentration or the 
coupling strength reduces T^.. This competition is influenced strongly by the number of 
states removed from near the Fermi surface by the formation of the superconducting state. 
In fact, Withoff and Fradkin 0] have pointed out that for the normal state analogue of 
this problem, in which the conduction band is assumed to behave as a semimetal, with 
density of states varying as A^(ci;) ^ |u;|'', r > 0, there is a critical coupling strength J^ 
below which the Kondo effect no longer occurs. One might expect the consequences of such 
a transition, were it to take place, e.g., as a function of temperature in an unconventional 
superconducting state providing the power-law density of conduction electron states, to be 
dramatic. In a previous study of this problem, however, we have shown that under most 
circumstances a finite density of Kondo impurities in the thermodynamic limit induces a 



finite density of states at tlie Fermi level, ensuring that a Kondo effect always occurs, but 
with reduced effective Kondo temperature. Nevertheless, unusual pairbreaking effects are 
still to be expected in unconventional superconductors due to the dynamical reduction of 
the Kondo scale caused by pair correlations. 

While this work was begun in anticipation of applications to heavy fermion supercon- 
ductors, recent measurements of penetration depth, |^ photoemission, [^ and Josephson 
tunneling on YBCO and BSSCO have provided evidence that the copper oxide super- 
conductors may be unconventional, possibly d-wave as well. This conclusion remains con- 
troversial, however, and experimental tests to distinguish conventional from unconventional 
pairing are of great current interest. One would like, for example, to develop an under- 
standing of the effect of doping Zn and Ni impurities in the Cu02 planes. Simple models 
of Zn and Ni acting as strong and weak potential scatterers in a d-wave superconductor, 
respectively, are consistent with some experiments at low doping levels, [|,B|,|10| but inconsis- 



tent with other measurements. |11,^2[ Since in some cases Zn impurities appear to possess a 



magnetic moment at higher temperatures, |]T2[ it is of interest to explore whether an s-d type 



exchange coupling of conduction electrons to an impurity embedded in an unconventional 
superconductor can describe the range of behavior observed. 

In this work we focus on basic thermodynamic and transport properties of supercon- 
ductors doped with Kondo impurities. Our aim is to develop a tractable, self-consistent 
scheme for the calculation of all basic properties of superconductors using methods known 
to successfully describe the most difficult aspect of the problem, namely the dynamics of 
the Kondo impurity. For this reason we have adopted the large- A^ "slave boson" approach 
of Barnes, [^ Coleman, H^ and Read and Newns, [0 as this approach is well-known to 



provide a good description of the spectral properties of the Kondo impurity for sufficiently 
low temperatures. We emphasize that this description is adequate for most of our purposes 

{T(.<Tk) even in the case of spin degeneracy N = 2, since the impurity spectral resonance 
is located at the proper position, i.e. exactly at the Fermi level in this approximation. 



II. FORMALISM 

We use the large- A^ slave boson technique for the SU{N) Anderson model describing an 
A^-fold degenerate band of conduction electrons, Ckm, m = 1, ...N with energy ek hybridizing 
through matrix element V with a localized impurity state fm- In general this Hamiltonian 
contains a term with the Coulomb repulsion U between two electrons present at the impurity 
simultaneously. In many compounds U is large (5-10 eV in the lanthanides) and for the 
purpose of studying the low temperature physics we will assume U = oo. The on-site 
repulsion term is then absent, but a constraint is added to ensure that the system remains 
in the physical part of the Hilbert space. The conduction band is assumed to have a constant 
density of states in its normal state, N{uj) = 1/2D = Nq. We also include a BCS-like pairing 
term of electrons on opposite sides of the Fermi sphere, 

H =Y1 ^kC^kmCkm + EfJ^ Pmfm + VY^[c^kmfmb + h.c] 
k,m "m k,m 

+ ^[A(A;)ctfc„ct_fc_^ + h.c] + A(5] f ^f^ + h^h), (2.1) 

k,m rn 

where A is a Lagrange multiplier enforcing the constraint nj + ri;, = X^m P mfm + i>^b = 1, 
preventing double occupancy of the impurity site. In the limit Ef —>■ — oo, NqV^/Ej =const, 
Eq. (2.1) reduces to the Coqblin-Schrieffer Hamiltonian with pairing studied in Ref. 0. 
Here we have chosen the more general form (2.1) to study deviations from single occupancy, 
Uf ^ 1, although we do not attempt to explore the fully developed mixed valent regime. 

The mean field approximation to this model, with mean-field amplitude (b) , leads to the 
two equations 
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dujf{co)-Ti (ro + r3)(G°(A;,cu + zO+)G/(cu + zO+)) , (2.3) 
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which determine (b) and e/, the latter being the position of the resonant state. Eqs. (2.2) 
and (2.3) should be solved self-consistently together with the gap equation. 



A(A;) = -TY.Y.Vkk'Tr^G{k',u^). (2.4) 

l^n k' ^ 

The full conduction electron Green's function G and the impurity Green's function G/ 

IGF\ [Gf Ff\ 

G = and G/ = (2.5) 

VFt G*J \Ff^ Gf*J 

are calculated from the diagrams shown in Figure 1, yielding 

G(cc;)-' = G\ujy' - S(cu) = cDro - e^ra - A{k)n, (2.6) 



Gfiiuy^ = Gj(cj) ^ - S/(cu) = euro - e/rs - An. (2.7) 

The Green's functions are now averaged over impurity positions in the usual way. The 
renormalized frequencies are calculated self-consistently from the Dyson equations, 

u = uj + au!/{-ui^ + e^f + A'^), (2.8) 

and 

u; = uj + T/i:o/[A^{k)-Cj^y^^y. (2.9) 

Here a = NnTTco/^n , n = h/NqTco is the scaled impurity concentration, and (...)^ is a 
Fermi surface average. The off-diagonal renormalizations are 

A(A;) = A(A;) + aA/{-LO^ + e^j + A^), (2.10) 

and 



A = r(A{k)/[A\k)-u^y ). (2.11) 

In superconductors with order parameters where the Fermi surface average in Eq. (8) 
vanishes, off-diagonal corrections vanish and A (A;) = A{k). This class includes but is not 
limited to odd-parity superconducting states. 

The energy scale F is the renormalized resonance width, F = {b)7rNQV^. The low- 



temperature Kondo scale in the large- A^ slave boson theory is given by Tk = a/F^ + e^ 



Although the width of the actual spectral feature corresponding to the Abrikosov-Suhl res- 
onance is modified below Tc, in what follows we will normalize all quantities with respect 
to this Tk, evaluated from Eqs. (2.2) and (2.3) with A = at T = 0. In the regime of 
principal interest, Tk > Tc, corrections to this definition will be small in any case. 

III. RESULTS FOR AN S-WAVE SUPERCONDUCTOR 

The early history of the problem of a Kondo impurity in an s-wave superconductor 
has been reviewed by Miiller-Hartmann. [|l^ Abrikosov and Gor'kov first discussed the 
pair-breaking effects of magnetic impurities weakly coupled via exchange interactions to 
conduction electrons. [Q Shiba |T^ extended this approach to treat strong scattering by 



classical spins, using the t-matrix approximation, and showed the existence of bound states 
in the gap. At finite concentration of impurities the bound states were found to form an 
impurity band whose width and center scaled with impurity concentration and exchange 
strength. These early works neglected the dynamical screening of the localized spin by the 
conduction electron gas. These effects were incorporated by Miiller-Hartmann and Zittartz, 
| 18| adopting an equation of motion decoupling scheme previously used by Nagaoka ll^ 



to calculate the dynamical spin correlations in the normal state. This approach correctly 
reproduced results in the Abrikosov- Gor'kov limit, Tk/T^. — > 0, and made the remarkable 
prediction of a reentrant superconducting phase if Tk -C Tc, subsequently observed in ex- 
periments on Lai_xCexAl2. pO| , ^ The failure of the decoupling scheme used to capture the 
correct crossover to Fermi liquid behavior in the normal state as T ^ invalidated the 
Miiller-Hartmann-Zittartz approach in the low temperature regime Tk > Tc, however. 

The physics of the Fermi liquid regime, Tk/Tc — > oo, was studied by Matsuura, Ichinose 
and Nagaoka |^ and by Sakurai [^ by extending the Yamada-Yosida theory ||2^ to the su- 
perconducting state. They obtained an exponential Tc-suppresssion with increasing impurity 
concentration n, Tc ~ Tcoexp{—pn/X), where A is the BCS dimensionless coupling constant, 
and p is a constant of order unity. This is commonly referred to as "pair- weakening" as op- 
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posed to pair-breaking, since the effective superconducting coupling constant is reduced due 
to correlations on the impurity site. The exponential form breaks down for concentrations 
sufficiently close to a critical tt-c, for which Tc = 0. In this regime the reentrant behavior 
found by Miiller-Hartmann and Zittartz does not occur. A further characteristic signature of 
the Fermi hquid regime is the reduced specific heat jump, C* = {AC / ACq) / (Tc/Tco))\Tc=Tco 
which is always less than one, [^,^ in contrast to the high temperature regime. 

Not surprisingly, qualitatively similar results were obtained by other early workers for 
Kondo and Anderson impurities using a variety of other approaches. [|l^] More recent treat- 
ments include the use of a self-consistent large- A^ p7|,|3[| , Monte Carlo p8[ and NRG methods. 



Schlottman |3^ has treated the mixed-valence regime using Brillouin-Wigner perturba- 
tion theory. Out of these efforts has evolved a qualitatively consistent picture of the effect 
of Kondo impurities on the superconducting transition, [^ but little understanding of the 



low-temperature properties of Kondo-doped superconductors because of the difficulty of the 
calculations involved. In this section we show that the current theory reproduces the known 
effects of Kondo impurities on the critical temperature, specific heat jump and bound states 
spectrum of an s-wave superconductor. 

A. Critical Temperature 

The simplest and most direct effect of impurity scattering on a superconductor is the 
suppression of the critical temperature. Scattering from impurities with internal quantum- 
mechanical degrees of freedom leads to deviations from the classic Abrikosov-Gor'kov predic- 
tion for the dependence of T^ on impurity concentration. [|l| These effects depend sensitively 
on the low-energy behavior of the self-energy S(c<j), which enters the linearized gap equation, 
obtained from Eq. (2.4) near Tc, 

ln(TJTco) = 27rT, V \ , ,, - V ^-—, (3.1) 

where B{uJn) = (tOn + r)^ + e/^ . 



The slope of the Tc-suppression evaluated at n = is 

[t,o dn ^^„ 27r ^^ {n + l/2)i?K) " ^ " ^ 

In an insert to Figure 2, we have plotted a numerical evaluation of Eq. (3.2) for an s-wave 
superconductor. Note the curve is drawn with a broken line for small Tk/Tco to reflect 
the fact that the slave boson mean field theory breaks down there. The maximum of the 



Tc-suppression is found to occur at Tk — 3Tc, similar to the NCA result. |23| The early 



high-temperature theory of Miiller-Hartmann and Zittartz |]18[ locates this maximum at 



Tk/Tc — 12, whereas in the more recent Monte Carlo calculation p8[ for the symmetric 
Anderson model with finite U the maximum slope of the Tc-suppression is at Tk — Tc. 
Unfortunately a direct quantitative comparison with the latter work is not possible, as the 
simulation is not performed in the fully developed Kondo regime. 

The present theory predicts an exponential decrease of Tc at small concentrations, in 



agreement with other theories of the Fermi liquid regime, p5| , p6 
Tc ^ Tcoexp 



" vl>(T;,/27rTco) 



Tk' 



c^TcIi-^HTk/Tco)], (3.3) 



K 



where \l/ is the digamma function. In Refs. [Q and |^3[ the initial suppression of Tc is 
proportional to ln^(Tx/Tco). A full evaluation of Eq. (3.2) for arbitrary concentrations and 
various values of the ratio Tk/Tc is shown in Figure 2. It is interesting to note that the 



theory reproduces the reentrant behavior characteristic of the high temperature regime |[T6| 
although we do not expect the theory to be accurate in this case (dashed line). 

As is evident in the Figure the current theory predicts no critical concentration Uc for 
which Tc = 0. This is a subtle point discussed by Sakurai, ||2^ who suggests that a failure 



to include the dynamics of magnetic scattering by states close to the Fermi surface can 
lead to such an effect. Such processes are included in the finite-[/ perturbation theory 
through Coulomb vertex corrections to the impurity averaged pair correlation function. In 
our U = oo theory, such vertex corrections arise first in leading order 1/A^ corrections due 
to the exchange of slave bosons, whose dynamics are neglected in this work. We expect that 
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effects arising from the absence of these fluctuations in the theory will be quantitatively 
small for T^^Tk, except for impurity concentrations so large such that Tc ^ T^q. 

B. Bound states 

For conventional superconductors, A{k) = A, Eqs. (2.9) and (2.11) are simply 

^ = ^ + r .. "" .,/, , (3.4) 

and 



A = r ^, ^ ,,, . (3.5) 

It is obvious from Eq. (3.5) that a finite gap in the conduction electron spectrum induces 
a gap in the impurity spectrum. When the Abrikosov-Suhl resonance, which develops at 
temperatures below T^, falls in the superconducting gap, bound states appear in the con- 
duction electron spectrum. These peaks are placed symmetrically relative to the center of 
the gap and their spectral weight depends on Tk/Tc and on impurity concentration. The 
density of states for N = 2 and A^ = 4 is shown in Figure 3. In the high-temperature regime 
the bound states emerge from the edges of the gap and move towards the center of the gap 
as Tk/Tc increases. [^ In the low-temperature hmit, the bound states disappear into the 
gap edges again. [Q Recently, Shiba et al. [^ studied the position of bound states in the 
gap in an s-wave superconductor using the numerical renormalization group (NRG). Their 
results cover both the high- and low-temperature limits reliably. The Monte Carlo study by 
Jarrell et al. [Q confirms the overall dependence of ub on Tk/Tc- 

Bound states correspond to the poles of the t-matrix for conduction electrons, T(u;) = 
V'^Gf{uj) and are given by 

, , /^ 2 r2a;2 y/2 r|^| , , 

\uj\ = \Tk + -r^ ^ TTT- (3.6) 

We compare the solution of Eq. (3.6) with the NRG result in Figure 4. There is a good 
agreement for Tk > Tc. For Tk ^ Tc the position of bound states ujb/A is a quadratic 



function of Tk/Tc, \ujb/ ^\ — 1 — 2A^/T^i^. The discrepancy between the NRG result and 
our calculation for Tr- < T^. is not surprising since the slave-boson theory fails in the high 
temperature regime. The slave-boson mean-field amplitude vanishes around T ~ Tk and 
the theory is unable to describe the crossover to temperatures above Tk- In Figure 4 we also 
show the spectral weight of the bound states in the gap which again agree with the NRG 
calculation for Tk ^ Tc- |^ 

To calculate the specific heat jump we use equations (2.5), (2.7) and (3.1), (3.2), and 
expand the gap equation near Tc in terms of A/T, 

TcQ _ _ 

n>0 V"--r -L/^A-^ -r"/-^V^n;; n>0 



-T.TTTTTj—mT-V. + 'Lzr^.A'^ril^t (3-7) 



'° T h (n + 1/2)(1 + a/B(w„)) ^ ,^„ « + 1/2 ^ 2°' birrj 



where 



and E{uJn) = 1 + a/ B{ujn). After expanding the free energy to fourth order in A we obtain 

c.(r.)-c„ra^?^(i-E^i%±^)'. (3.9) 



As can be seen in Figure 5, the dependence of AC/ACq on T^/Tco, for Tk^Tcq, is qual- 
itatively different from the Fermi liquid regime. In that limit, C* c^ 1 — l/ln(Tft:/2Tco), 
see Figure 6. Ichinose [^ and Sakurai |^ have obtained a somewhat different result. 



C* ~ 1 - l/ln^(r;,/T,o). 

IV. IMPURITIES IN UNCONVENTIONAL SUPERCONDUCTORS 

The problem of an Anderson impurity in an unconventional superconductor is of con- 
siderable interest for several reasons. First, it serves as a testing ground for the ideas of 
Withoff and Fradkin, who argued that in the analogous Kondo problem with power law 
conduction electron density of states, N^u) = C\uj\'^, there existed a critical coupling below 
which impurities are effectively decoupled from the conduction band. BJ36[| We showed in 
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Ref. that in an unconventional superconductor a critical coupling indeed exists for the 1- 
impurity problem, but that for a finite density of impurities there is always a finite density of 
conduction electron states at the Fermi level, provided A^ = 2 or if r < 1. This is consistent 
with phenomenological studies of potential scattering in unconventional superconductors. 
P7| , |55[| This does not exclude the possibility of observing some vestige of this transition as 
the Kondo temperature becomes quite small, however. 

Secondly, Kondo impurities in unconventional superconductors have been proposed as 
analogues of defects in Kondo lattices, [^ with the argument that a vacancy in a Kondo 



lattice may induce a relative phase shift close to 7r/2. While such a picture clearly cannot 
account for the degradation of coherent transport by such defects, it is interesting to further 
explore the analogy, which has been successful in describing the superconducting state of 
heavy fermion materials. 

Finally, there is currently considerable interest in the possibility that the cuprate su- 
perconductors may be characterized by an unconventional order parameter, as suggested 
by both experimental [§|,§|,0!|§| and theoretical [^ studies. Studies of YBa2(Cui_xMx)307, 



where M are Zn or Ni dopants in the CuO planes, suggest that the behavior of the two ions 
can be quite different. Maple and collaborators have claimed that Pr may act as a Kondo- 
like dopant in YBCO. We believe that a Kondo- or Anderson-type model may be sufficiently 
rich to describe the unusual features of the Zn and Ni doping studies as well, in conjunction 
with an understanding of the different gap nodes in unconventional superconductors. To 
this end we study the effects of Kondo impurities on model unconventional states, bearing 
in mind that the theory may be inadequate to describe the fully developed magnetic regime 

We chose to do our calculations for the axial, A(fc) = Ao{kx + iky), and polar, A{k) = 
Aq/Cz, states for simplicity. These two states, nominally p-wave pairing states over a spherical 
Fermi surface in 3D, are quite generally representative of two classes of order parameters. 
States with order parameters vanishing at points on the Fermi surface, like the axial state. 
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have low temperature properties associated with a quasiparticle density of states N^u) ~ tu^, 
whereas the states with hnes of nodes, hke the polar state, correspond to A^(u;) ~ uj. More 
complicated order parameters, such as those having a rf-wave symmetry, will have similar 
properties at low temperatures, since the main factor determining the low-T properties is 
the order parameter topology, i.e. whether there are points or lines of zeros of the order 
parameter. In the case of points (lines) of nodes, the low-temperature specific heat of pure 
superconductors is proportional to T^ (T). The deviation of the penetration depth AA from 
its zero temperature value A(0) along the main axes of symmetry of the order parameter 



is either ~ T or ~ T in the axial state, according to the direction of current flow, ^ 
whereas for the polar state it is either ~ T or ~ T^. The presence of strong impurity 
scattering complicates the picture and these power laws do not hold in general. 

Due to the absence of the Anderson theorem for p-wave-like or d-wave-like superconduc- 
tors, the influence of impurity scattering on the critical temperature is qualitatively different 
from that of s- wave-like superconductors. P2|J^ For the unconventional states of interest, 
there are no off-diagonal corrections to the superconducting Green's function, A = A. In 
this case T^ is determined by 

MTjTco) = 2tiTc y — -, , , "^,, — - - y ^. (4.1) 

The initial Tc- suppression is then given by 

and approaches —NT'^/8T^, in the limit Tk/Tco — > cxd, see Figure 7. Note that cJn/A -^ 
aT/AT^ , as Tc —> , and as a consequence the critical concentration ric is finite. 

The specific heat jump at T^. can be found by the method already mentioned in the 
preceding section, 

Svr^iVoTc 



Csin) - Cn{T,) 



h 
where 
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>r- / G{un) 

„Vo V((" + V2)i5K) + «/5(^n))^ 



(4.3) 



with b = 4/5 (3/5), and c = 2/3 (1/3), for the axial (polar) state, and H{iJn) = 
aT [2{iJn + ry/Biun) - 1] /Bioun) , GM = 2auOn{oOn + T){n + 1/2 + T/2TxT,)/B\uOn) + 
ar/27rr,fi(cu„) . 

The effect of Kondo impurities on the density of states is shown in Figure 8. As in 
the s-wave case, the resonant states move towards the edges of the gap as Tk/Tc increases, 
except in the special case N = 2, where they are pinned at the gap center, f^ In all cases the 
impurity bands are broadened relative to the s-wave case by the continuum in which they 
are embedded. In the limit Tk —>■ oo, the current theory coincides with the results given by 



phenomenological t-matrix treatments [p7||38| for N = 2. 



To obtain the specific heat we differentiate the entropy with respect to temperature 
C = TdS/dT. The entropy is given by 

POO 

S = -kB du;N{u;)[f\nf + {l-f)\n{l-f)], (4.5) 

Jo 

where / = f{uj) is the Fermi function. Note that the density of states N{uj) is calculated 
self-consistently, using Eqs. (2.2-2.4). 

The presence of resonances at low energies leads to pronounced features in the low-T 
specific heat, as shown in Figure 9. For T^ ^ T^, A^ = 2, these features are identical to 
those predicted by the phenomenological theory of Refs. [^ and ||38|. For smaller Tk/Tc, 



the resonance sharpens, as is evident from, e.g.. Figure 8. In this case resonances in the 
low-temperature specific heat may be quite dramatic. Figure 10 shows the temperature 
dependence of C/T for A^ = 2 in the axial state for the case when the resonance is very 
narrow, T <^ T^ and just above the Fermi surface, ej <^ T^.. Such pronounced features 
are possible only when the bare impurity level is close to the Fermi surface, and the hy- 
bridization is weak. They will be sharper for superconducting states with larger exponent 
r in the unperturbed density of states N{uj) ~ C'l^uT ^^ lo* ^- These anomalies may be 
observed experimentally at sufficiently low temperatures. In this context, it is interesting 
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to note that a sharp peak in C /T has been observed in the heavy fermion superconductor 
UPt^ at 18 niK. [^ This peak is present at roughly the same position also in the normal 
state at magnetic fields B > Bc2, as might be expected in a situation where the Kondo 
temperature is significantly smaller than the critical temperature. If such an interpretation 
of the measurement of Schuberth et al. in these terms is correct, we would expect the size 
of the peak to scale with other measures of the defect concentration, such as Tc the size of 
the specific heat jump, in different samples. We note, however, that fields of order 1 Tesla 
would normally destroy a many-body resonance of the usual magnetic type. 

We now discuss the low-temperature response functions of the superconductor in the 
presence of Kondo impurities, which have not to our knowledge been previously calculated 
in the strongly interacting Fermi liquid regime of interest. The effect of Kondo impurities 
on the electromagnetic response has been calculated in the phenomenological t-matrix ap- 
proach mentioned above, and used to analyze experiments on heavy fermion and high-Tc 
superconductors with impurities, but no microscopic theory is available. The London pene- 
tration depth is obtained from A = ——7^(0, 0) , where K(0, 0) is the electromagnetic 

response kernel. The kernel is given by the linear response formula 

ir*^(0,a^) = -- / deY,Tr{hk,Gik,Un)Gik,LU,,-nm))k. (4.6) 

ZTYIC J —oo ji 

In the static limit, Qm ~^ , the kernel becomes 

A-«(0. 0) = -5^ E / ^^.*,.^!W . (4.7) 

mc '^ J 4:7!- (A^(/i;) + tu^)'^/^ 

where the integral represents the angular average. Here we specialize again to the N = 2 
case. The temperature dependence of (A^^)^^ in the polar state, which corresponds to the 
component of the superfluid density within the plane containing the line of zeroes of the order 
parameter, is shown in Figure 11. The low-T behavior changes from linear to quadratic upon 



doping. A similar result was obtained earlier PT| , ^ within the phenomenological theories. 
The results for T^ ~ Tc are qualitatively very similar to those shown. Finally, A^^ at 
T = along the main axes of symmetry for the axial and the polar state as a function 
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of impurity concentration is presented in Figure 12. We note that the largest component 
of the penetration depth scales in the case of line nodes as X{n) / Xl{0) — 1 ~ n^/^ at low 
concentrations, and A~^(n) ~ log{nc/n) close to the critical concentration. It is clear that 
both the concentration and temperature dependence of the penetration depth components 
may be important tests of gap anisotropy. 

V. CONCLUSIONS 

We have presented a slave boson theory of Kondo impurities in superconductors which has 
the advantages of being applicable in the Fermi liquid regime and being relatively easy to use 
in calculating quantities of experimental interest at all temperatures in the superconducting 
state. The theory has been shown to reproduce all of the qualitative features of the physics 
found by previous theories for large T^, and is asymptotically in quantitative agreement with 
"exact" NRG calculations for the s-wave case. It is furthermore capable of going considerably 
beyond currently available theories in that practical calculations of superconducting response 
functions are possible at low temperatures in the superconducting state, and can be easily 
generalized to the unconventional states of great current interest in the heavy fermion and 
high-temperature superconductivity problems. 

The single exception to this success is the failure to properly describe of the s-wave 
superconductor at very large impurity concentrations, in that the theory predicts an infinite 
critical concentration. In practical terms this is quite academic, as all independent impurity 
analyses will break down due to interimpurity interaction effects long before any putative 
critical concentration is achieved. Nevertheless, this formal shortcoming of the theory exists 
and must be addressed. Preliminary analysis of fluctuations about the saddle point has 
convinced us that the Gaussian fluctuations to the scattering amplitude arising from the 
slave boson dynamics will be sufficient to induce a critical concentration, in analogy to the 
works of Matsuura et al. [^ and Sakurai [^, and that the theory as it stands is sufficient 
to describe the Fermi liquid regime everywhere except for very large concentrations in the 
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s-wave case. The difficulty does not arise in the unconventional case. 

In unconventional superconductors, we have shown that the phenomenological theories 



of Refs. [0 and ||38| will be reproduced by the microscopic theory presented here in the 
limit N = 2 and Tk/Tc —>■ oo. The decondensation of the slave boson amplitude prevents 
extension of the theory into the high temperature regime, but qualitatively it is clear that 
the effect of lowering the Kondo temperature is to sharpen the many-body resonance near 
the Fermi surface, but lower its weight. Thus the effect of resonant scattering leading to 
low-energy gapless effects in superconducting thermodynamic and transport properties is 
reduced. Impurities with larger orbital degeneracy may lead to similar resonances away 



from the Fermi level, possibly similar to those observed by Maple et al. |^ in specific heat 
experiments on Pr-doped YBCO. 

The effect of Zn and Ni-doping on superconducting properties of YBCO suggests that 
these two dopants may be described within a traditional Kondo-type picture in conjunction 
with the theory discussed here. Ni, with spin one, may possibly be treated as a higher- 
degeneracy scatterer below its Kondo temperature. As we have seen, such an impurity acts 
as a weak scatterer compared to the N=2, T^ ^ T^ resonant scattering case, a possible 
model for Zn. On the other hand, recent NMR measurements appear to suggest that a 
moment forms around the Zn site. Within an isolated impurity picture, this would suggest 
that Tx <^ Tc, where the present theory is not applicable. A more plausible explanation, 
however, is that the local spin correlations induced by a missing Cu must be accounted 



for, as suggested by Poilblanc et al. [Q Further experimental and theoretical work on this 
problem is clearly essential. 
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FIGURES 
FIG. 1. Dyson equations for the conduction electron and impurity Green's functions. 

FIG. 2. The critical temperature for an s-wave superconductor as a function of impurity con- 
centration. The inset shows the slope of this dependence evaluated at Tc = Tcq. 

FIG. 3. Conduction electron and impurity spectral functions in the Kondo limit in an s-wave 
superconductor for N = 2 and A^ = 4. The solid and dashed lines correspond to Tr- = T^. and 
Tft- = lOTc, respectively. The concentration of impurities in all cases is n = 0.4. 

FIG. 4. The position and spectral weight of the bound states in the gap for an s-wave supercon- 
ductor with Kondo impurities. The solid (dashed) line is the location (spectral weight) of bound 
states. Only one of the bound states is indicated here, the other one is located at positive energies, 
symmetrically with respect to the gap center. Circles (triangles) refer to the position (spectral 
weight) obtained from an NRG calculation. [29 

FIG. 5. Specific heat jump as a function of Tc/Tcq. We included the result of the 
Abrikosov-Gor'kov theory [45| for comparison. 



FIG. 6. The derivative of the specific heat jump evaluated at zero impurity concentration as 
a function of T^/Tc for an s-wave superconductor. The dash-dotted line shows the asymptotic 
behavior, 1 — l/ln(Tft'/2Tco). The dashed line is the asymptotic form at Tk/Tc — > oo found in 
Refs. |2§ and H, C* ~ 1 - l/\n^{TK/Tco). 



FIG. 7. Critical temperature vs. impurity concentration for unconventional states considered 
in this work. The initial slope at Tco is shown in the inset. 

FIG. 8. Conduction electron and impurity spectral functions in the Kondo limit in a polar 
state for N = 2 and A^ = 4 . The solid and dashed lines correspond to Tk = Tc and Tk = lOTc, 
respectively. The concentration of impurities is n = 0.4 
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FIG. 9. The low-T part of C /T in the polar state for A^ = 2 in the Kondo hmit. The inset 
shows C /T for n = O.Olnc over the full temperature range. 

FIG. 10. The low-T part of C /T in the axial state. Note the sharp resonance in the density of 
states at low energies (the inset). 

FIG. 11. The inverse square of the penetration depth vs. temperature at several concentrations 
of impurities in the polar state. 

FIG. 12. The inverse square of the T = penetration depth in the two principal directions as 
a function of impurity concentration for the polar (full lines) and the axial state (dashed lines). 
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